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Abstract 



O ' We further discuss possible sgoldstino interpretation of the obser- 

vation, reported by the HyperCP collaboration, of three S + — ► pfi + [i~ 
Qh| decay events with dimuon invariant mass 214.3 MeV within detector 

O ■ resolution. With sgoldstino mass equal to 214.3 MeV, this interpreta- 

O , tion can be verified at existing and future B- and ^-factories. We find 

that the most natural values of the branching ratios of two-body B- 
and D-meson decays to sgoldstino P and vector meson Y are about 
^ . 10~ 6 -T- 10~ 7 . The branching ratios of 0-meson decay <j> — > Pj is es- 

timated to be in the range 1.8 • 10 -13 -j- 1.6 • 10~ 7 , depending on the 
hierarchy of supersymmetry breaking soft terms. Similar branchings 
for p- and cj-mesons are in the range 10 -14 -j- 3.4 • 10~ 7 . 

Recent publication [1] of anomalous events 

E — >pfi + fi~ , (1) 

observed in the HyperCP experiment, renewed interest to models with light 
particles whose interactions with the Standard Model fermions violate flavor. 
One of the well-motivated classes of models of this kind, supersymmetric 
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models with light sgoldstino 1 , has been suggested by the HyperCP collabo- 
ration [1, 9] to explain the anomalous events. The latter are three events (1) 
where all dimuon masses fall into a single bin within the detector resolution; 
notably, no more decays (1) have been observed. This coincidence has been 
suggested as the first evidence for sgoldstino production in two-body decay 
with subsequent decay of sgoldstino X into yU + /x~-pair, 

E — > pX , X — ► • (2) 

Sgoldstino mass has been fixed from measured dimuon masses as 

m x = 214.3 ± 0.5 MeV . (3) 

It was pointed out in Refs. [10, 11, 12, 13] that, generally, X-particle can 
be either pseudoscalar or pseudovector, but neither scalar nor vector. Thus, 
only pseudoscalar sgoldstino (denoted by P in what follows) in models with 
parity conserving sgoldstino interactions (see Ref. [14] for a description) is vi- 
able [12]: light scalar sgoldstino, as well as pseudoscalar sgoldstino in models 
with parity violating sgoldstino interactions, cannot explain the anomalous 
events because of the strong upper bounds on the widths of two-body kaon 
decays, K — > 7iX(X — > 

From the presented [1] branching ratio of the decay (1) (which is some- 
what above the Standard Model predictions [15, 16]), the following product 
has been estimated [12], 

\h$\ ■ BtV 2 (P -> /iV) = 3.8 • 10- 10 . (4) 

Here is sgoldstino coupling constant to pseudoscalar current d^s, 

C Pds = -P ■ (M? -di^s + h.c.) , 

and Br(P — > denotes sgoldstino branching ratio to muons. 

Additional information came from the geometry of the detector and the 
measured energies of muons, which enables one to put an upper limit on 
sgoldstino lifetime [12], 

r x < 2.5 • 10~ n s . (5) 

1 See, e.g., Refs. [2, 3, 4, 5] for examples of models and Refs. [6, 7, 8] for the description 
of sgoldstino interactions. 
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Note that a lower limit [12] on sgoldstino lifetime coming from the current 
limit on supersymmetry breaking scale is about 



r x >1.7-l(n 15 s (6) 

for sgoldstino mass (3). 

The suggested explanation of HyperCP anomalous events awaits a test in 
other experiments. The processes which can be used to definitely confirm or 
rule out the sgoldstino explanation have been discussed in Ref. [12]. These are 
three-body kaon decays 2 and direct sgoldstino production in e + e~ collisions. 
The rates of these processes are governed by the same parameters of the 
sgoldstino Lagrangian, which are constrained by Eqs. (3), (4), (5) and (6). 

The parameters are not directly constrained by the HyperCP result. 
These are other flavor- violating sgoldstino coupling constants, hji (h^\ j, I = 

1,2,3, and hff\ except for h^), entering 

C P = -P ■ (h jt ■ ] 3 i^U + h.c.) , (7) 

where fj are fermions of the Standard Model, and sgoldstino flavor-blind 
couplings which do not contribute to sgoldstino width because of kinematical 
constraints imposed by the small sgoldstino mass (3). Nevertheless, these 
couplings make it possible to extend searches for sgoldstino. Though the rates 
of corresponding processes cannot be predicted because of lacking knowledge 
of model parameters, one can estimate the most natural ranges for these 
rates by making use of the expected patterns of MSSM soft terms. With 
sgoldstino mass being fixed, these predictions can be tested experimentally, 
giving an opportunity to find additional evidence for light sgoldstino (but, 
generally, not to compromise the proposed sgoldstino explanation of the 
HyperCP result). This issue is discussed in this letter. 

The dominant 3 sgoldstino decay modes are [14] P — > ji + ji~ and P — > 77. 
A natural range of values of the ratio of their widths is 



2 K — ► irirX, K — > and f2~ — > S~ X decays have been discussed in Refs. [10, 

11, 12, 13] as signatures to test the hypothesis (2) with some pseudoscalar or axial vector 
particle X in a model-independent way. 

3 The kinematically allowed decay mode P — > e+e~ is typically suppressed if sgoldstino 
mass is about 214 MeV. 
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However, much larger values for this ratio, such as 10 4 are not excluded (see 
[12] for discussion). So, we are interested in rare decays with fi + fi~ or 77 
with invariant mass 214 MeV in final states. 

As promising processes, we consider two-body B- and D-meson decays 
(to probe flavor-violating sector of sgoldstino couplings) and two-body light 
neutral vector meson decays (to probe flavor-blind sector of sgoldstino cou- 
plings). These searches may be performed at B-factories where they would 
cover the most natural part of the relevant parameter space, while searches 
at ^-factories seem less promising, as one can probe there only a part of the 
relevant parameter space leaving the rest for future experiments. 

We begin with sgoldstino flavor-violating interactions (7) in quark sector. 
The coupling constants hij entering these interaction terms are determined 
by the ratios of the elements of left-right part of squark squared mass matrix 
rnif^ 2 to vev of auxiliary field in goldstino supermultiplet F, the latter being 
of the order of the squared scale of supersymmetry breaking in the complete 
theory, = m^ R ^ 2 j '(V2F) . In particular, non-zero off-diagonal entries of 
squark squared mass matrix lead to two-body pseudoscalar B- and D-meson 
decays to a vector meson V and pseudoscalar sgoldstino 4 , 

Pb,d -> VP ■ (9) 

Let p, p~y be momenta of Pb,d and V, respectively, q = p — p r , and ey 
be the polarization vector of V. With these notations, the hadronic matrix 
elements 

°P B ,D,y (<? 2 ) = (V{pv,e^)\qjl 5 qi\PB,D{p)) 

entering the decay amplitudes M (Pb,d — > VP) = hj t ■ P (q) ■ Pb Dt y (q 2 ) 
(where P(q) is the wave function of the outgoing sgoldstino) can be ex- 
pressed via experimentally measured and/or theoretically predicted form- 
factors A^ B ' D,r \q 2 ) as follows (see, e.g., Ref. [17] for notations and other 
details), 

PB . D .r (f) = 4 PB '°''*V) ■ , (10) 

rrij + mi 

where my and m^ are masses of V and quark q iy respectively; we use m^ = 
4.8 GeV, m c = 1.35 GeV, m s = 0.13 GeV [17] in our estimates. Making use 

4 Two-body decays of pseudoscalar mesons into pseudoscalar mesons and pseudoscalar 
sgoldstino are strongly suppressed by parity conservation in sgoldstino interactions. 
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of the expression (10), one finds the partial width 



hn\ 2 f ,(P*n.r), m p B ,n X3 ( m PB,Di m r,m P ) 



(11) 



and X(m u m 2: m 3 ) = ^ (l - {m2 ^ 3? ) (l - ( " 3r j ■ The values of the 

form factors A B ' D ' entering Eq. (11) can be found in literature [18, 19, 
21, 17]. 

To illustrate the level of precision required to probe the models consistent 
with HyperCP events, we present estimates of the widths in three different 
types of supersymmetric models. Two of them are 

model I : hji ~ , j ^ I 
A 

model II : hji ~ — • max{m^m\) , 
F 

where A is a flavor-independent constant. The third model is a concrete phe- 
nomenologically viable example of left-right supersymmetric model, where 
parity conservation in sgoldstino interactions is guaranteed [14]. In the first 
model all off-diagonal entries in squark squared mass matrix are of the same 
order. In the second model there is a hierarchy in matrix elements hji re- 
flecting the hierarchy of quark masses. The latter situation is more realistic, 
at it is typical for minimal supersymmetric models like mSUGRA or models 
with gauge mediation of supersymmetry breaking, where soft supersymme- 
try breaking trilinear terms are proportional to the Yukawa matrix. For the 
model III we take a general left-right SUSY model of the type presented in 
Ref. [14], with left-right symmetry broken at energy scale Mr = 3 TeV. We 
assume (see notations in Ref. [14]) universality at the SUSY breaking scale 
y/F = 30 TeV, i.e. A« = Y®A®, choose = A^ = M 1/2 , tan/3 = 3.6 
and neglect mixing between Higgs doublets in doublet- doublet splitting [20] . 
The left-right entries in the squark squared mass matrices rhj^J^ are ob- 
tained by making use of one-loop renormalization group equations for gauge, 
Yukawa, soft trilinear coupling constants and gaugino masses [22] . The value 
of A^ at M R = 3 TeV is fixed by (4). Under our assumptions, all relevant 
parameters in model III are then completely determined. 

For these three models it is straightforward to estimate the partial widths 
of the processes (9) with real sgoldstino decaying into fi + fi~, Pb,d —> yP(P — > 
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decay 


^7 


a(Pb n,r) 
^0 


Br ( mode i I) 


Br ( modc i II) 


Br (model III) 


B s - 0P(P - 
B s -> K*°P(P -> 


"■13 


0.42 [18] 
0.37 [18] 


6.5- 10- 9 
5.3- 10~ 9 


8.8- 10~ 6 
7.2 • 10~ 6 


8.7- 10- 6 
2.3- 10~ 7 


P+ -> D*+P(P -> 

5+ - j d:+p(p - 

P+ -> P*+P(P -> 


u 13 
h (D) 

h (U) 
n 12 


0.14 [19] 

0.14 a 
0.23 [21] 


3.2 • 10- 10 

3.0 • 10~ 10 

4.1 • 10~ 10 


4.4- 10~ 7 
4.0 • 10~ 7 
4.4- 10~ 8 


1.4- 10- 8 
4.0 • 10~ 7 
8.2 • 10~ 7 


P+ -> K* + P(P -> 
B° ^ K*°P(P ^ 


h (D) 
rl 23 


0.31 [17] 


3.8- 10~ 9 
3.5- 10~ 9 


5.2 • 10~ 6 
4.8 • 10~ 6 


5.1 • 10~ 6 
4.7- 10~ 6 


B° -> pP(P -> 
P+ -> p+P(P -> 


h (D) 
a 13 


0.28 [17] 


3.1 • 10~ 9 
3.3- 10~ 9 


4.2 • 10~ 6 
4.6- 10~ 6 


1.4- 10~ 7 
1.3 ■ 10 7 


P° -> pP(P -> 

p+ -> p+p(p -> /xV") 


h (U) 


0.64 [17] 


1.4 ■ 10- 9 
3.5- 10~ 9 


1.5- 10~ 7 
3.7- 10~ 7 


2.8 • 10~ 6 
7.0 • 10~ 6 





a We did not find any estimate of this formfactor in literature and use this value as an 
order-of-magnitude estimate, which is sufficient for our study. 



Table 1: Branching ratios of decays Pb,d — > yP(P — > in the models 

I, II and III. Branching ratios of decays Pb,d "VP{P — > 77) are given by 
the same numbers multiplied by T(P — > 77)/r(P — > 



. The results are summarized in Table 1. For the widths of similar pro- 
cesses, but with sgoldstino decaying into photons, Pb,d "VP{P — > 77) , 
one gets the same numbers multiplied by the ratio T(P — > 77)/T(P — > 
whose estimates are given in (8). 
Comparing the results presented in Table 1 with statistics of B- and D- 
meson decays collected by B-factories, one concludes that both fi + fi~ and 
77 decay channels can be probed for the most natural choice of parame- 
ters (models II, III). Moreover, a part of expected region for is already 
excluded by the results [23] 

Br (P+ -> rW) < 2.2 • 10- 6 , 

Br (B° -> K*°n + fi-) = (1.3 ± 0.4) • 10~ 6 . 1 ' 

Study of the model I requires, generally, higher statistics. At B-factories 
model I could be probed if sgoldstino decay mode into photons dominates by 
one to two orders of magnitude over mode. 

It is worth noting that coupling constants h[^\ and hfp determine 
the rates of several decays each. Hence, the ratios of the corresponding 
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rates do not depend on the values of these couplings. This would allow for 
independent check of sgoldstino interpretation, would any of the anomalous 
decays listed in Table 1 be observed. In particular, Eq. (12) implies that 
branching ratios of B s — > (f)P(P — > p + /i~) and £?+ — > D* + P(P — > p + p _ ) 
have to be at least three times smaller than the numbers in the two last 
columns of Table 1. 

Another signature of models with light sgoldstino we would like to discuss 
is neutral vector meson Y decays to sgoldstino and photon, 

r -> P7. (13) 

The most sensitive are the decays of 0-, u- and p-mesons, whose rates are 
almost saturated by sgoldstino couplings to gluons, as we find below. For 
J/ip- and T-mesons the corresponding branching ratios do not exceed 10~ 9 
(see Table 5 in Ref. [27]). 

7 

y 7 



a) 




Figure 1: Diagrams contributing to the decay V — > 7P. 

There are two different contributions, .Mi and M2, to the amplitude of 
the decay (13). These are due to sgoldstino couplings to photons and to 
gluons, 

A = 7^ AP {^PF^ + ^yPG^Gl)} , (14) 

(where M 77 = Mi cos 2 9 W + M 2 sin 2 6 l w and Mi, M 2 and M 3 are masses of 
C/(1) Y -, SU(2) w - and S , f/(3) c — gauginos, respectively) and sgoldstino flavor- 
conserving interactions with quarks (7). 

The relevant interactions of vector meson "V are V — 7 mixing [24], 

£ r i = -a f m 2 r i/, , ^ ~ 8.6 (15) 
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(where my is the mass of the corresponding vector meson and (a p , a U) a^) = 
(^/2, — |^), and coupling to photon and pseudoscalar mesons [25], 

Cy,2 = egy^e^xpd^dxyp^ , (16) 

where the values of coupling constants gy^ are tuned to saturate the 0? — > 
^7 or ~Y — > ^7 decay rates and are presented in Ref. [25]. 

Couplings (14) and (15) give rise to the following contribution to the 
decay amplitude 

Mi = e» • i-a y • ^ • f-iv^^ • e^p x k p el(k) , 
g m r \ t J 

where p and k stand for momenta of ^-meson and 7, respectively, while e(p) 
and e{k) denote their polarizations. This contribution is illustrated in Fig. la. 
The interactions (14) give rise to P — 8? mixing [27], 



where 



= C P » ■ P@> , (17) 



7cm n M 3 m u — m d 
3a s (M 3 )F m u + m d 

nM 3 m 2 ( l) Bf w Am s ( R . x 

vW) = V3a F vW) V6F V C ° S Sm ) ' 

fvW) = fa c °s #8 =F V2f sin 9 , 

where upper sign stands for 77- meson and lower for rj'. Here B w m 2 K /m s ~ 
2 GeV and we neglect small (~ 1CT 2 ) mixing between 7r° and 77, -//-mesons, 
77 — 77' mixing angle = —15.4° and we adopt two mixing angle scheme for 
the parametrization of the gluonic matrix elements [25], [28] 

( ^3q^ g ^ |0) = ^ m 2 w)fv(vi) (18) 



3a s (M 3 ) ~ m u -m d 2 

(7r I — ^ — ^c^|o> - -— T — /.m, . 

where we use /„ = 1.17/ w , / 8 = 1.26/* and # 8 = -21.2°, # = -9.2° [26]. 
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Couplings (16) and (17) give rise to the contribution 



M 2 = (~ W 7 ) (p) e^p x k p el(k) 



m p — mi 



■iC. 



Pr, 



illustrated by the diagram shown in Fig. lb. Here sum runs over 7r°, r] and 
i]' mesons. 

Finally, the rate of the decay Y — > P7 in the case of unpolarized vector 
meson reads 



r (r -> p 7 ) = 



r 2 



2 \i 
mp) 



96tt 



m: 



CV = 



77 



+ 



m. 



(19) 



In most models, the gluino mass is large, so that the dominant contribution 
comes from sgoldstino couplings to gluons. The corresponding coupling con- 
stant is proportional to the ratio M 3 /F. It reaches the maximal value in the 
unitarity limit M 3 ~ \[F at the smallest value of VF~, which is consistent with 
current bounds on the scale of supersymmetry breaking, \/F > 500 GeV. So, 
we have 

^ ~ 2 • 10- 3 GeV" 1 . (20) 

J max 

If M 3 /F is much smaller than this value, the dominant contribution to 
comes from either M 77 /F or A/F terms, at least one of which is bounded 
from below by the upper limit on sgoldstino lifetime (5), 



M 



77 



~ 1.5 • 10' 5 GeV~ 



(21) 



When considering the case of small M 3 , we present the numerical results for 
the case \A\ <C M 77 . In the opposite case one obtains about two orders of 
magnitude smaller numbers for the lower bounds of the intervals. 

Thus, the branching ratios of "V — > P7 decays are expected to be within 
the intervals 



1.8- 10 -13 < Br* < 1.6-10 



-7 



9.1 • 10~ iO < Br p < 3.3-10"' 
1.8 • 10~ 14 < Br ul < 3.4 • 10~ 7 



(22) 
(23) 
(24) 
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The produced sgoldstino subsequently decays either into photons or into 
muons. Comparing the values of the branching ratios in the interval (22) to 
the collected world statistics of 0-mesons (about a few billion) one concludes 
that the part of the interval (22) with the largest branchings may be tested 
with existing experimental data at ^-factories, while the large part could be 
tested only with significant increase in statistics. Note that observation of 
two-body decays of different vector mesons would offer an opportunity to 
measure all relevant parameters M 77 , A, M3, F. 

To conclude, by studying two-body decays of B-, D-, 0-, cu- and p- 
mesons one can probe the models with light sgoldstino, capable of explain- 
ing HyperCP anomalous events. Currently available statistics is sufficient 
to probe the most natural left-right symmetric models and models with 
the hierarchy in sgoldstino flavor-violating couplings similar to the hier- 
archy of quark masses. The rest of models awaits larger statistics. The 
same conclusion holds for quite similar baryon two-body decay modes like 
IT -> E~P [10, 11, 12], tt° c -> S°P, etc.: though the rates of (some of) these 
decays can be predicted from the HyperCP data, the current statistics is too 
low, so that future experiments are needed to test the sgoldstino explanation 
by searching for these decays. 

If observed, these decays would not only confirm the sgoldstino interpre- 
tation of HyperCP results (which have been claimed as the first evidence 
for supersymmetry [9]) but would give an opportunity to measure sgoldstino 
life-time and to probe the pattern of sgoldstino flavor violating couplings to 
quarks. The former would fix the ratios of MSSM soft terms to the squared 
scale of supersymmetry breaking, while the latter would enable one to esti- 
mate the values of off-diagonal entries in squarks squared mass matrix, which 
can be tested, in turn, by searching for FCNC processes. 

The negative results of searches for sgoldstino in two-body B- and D- 
meson decays at the level of branchings as high as 10~ 7 -^10~ 8 (models II and 
III) would imply either non-sgoldstino explanation of the HyperCP results 
or fairly special pattern of off-diagonal squark masses, atypical for simple 
super symmetric extensions of the Standard Model. Hence, some mechanism 
(e.g., additional flavor symmetry) should have to work in that case to provide 
for this pattern. 

We are indebted to A. Bondar, M. Danilov, S. Eidelman, P. Pakhlov, 
V. Rubakov for numerous helpful discussions. This work was supported 
in part by grant NS-7293.2006.2 (government contract 02.445.11.7370), by 
RFBR grant 05-02-1736 and by the Dynasty Foundation (awarded by the 
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